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I.   INTRODUCTION 

The  text  by  Hadley  and  Whitin  [Ref.  1]  is  a  standard 
reference  on  mathematical  inventory  models.   A  particular 
class  of  models  presented  therein  which  has  received  wide 
acceptance  is  known  as  the  lot-size,  reorder  point,  or 
(Q,r),  model,  where  r  is  defined  as  that  inventory  level  at 
which  a  reorder  for  a  quantity  Q  is  initiated.   Both  Q  and 
r  are  assumed  to  be  continuous  variables.   The  authors 
begin  their  treatment  of  stochastic  customer  demand  by 
deriving  two  models  referred  to  as  the  backorders  case  and 
the  lost  sales  case,  respectively.   In  the  backorders  case, 
all  demand  occurring  when  the  system  is  out  of  stock  is 
eventually  filled,  whereas  in  the  lost  sales  case,  this- 
demand  is  lost  forever. 

The  models  derived  by  Hadley  and  Whitin  consist  of 
long-run  average  annual  variable  cost  functions,  both 
labeled  K(Q,r).   Each  function  consists  of  a  sum  of  three 
cost  expressions,  the  first  representing  a  set-up  or  ordering 
cost,  the  second  an  inventory  carrying  cost,  and  the  third 
a  backorders  or  lost  sales  cost  as  appropriate.   The  authors 
assert  that  these  functions  are  convex  in  the  variables  Q 
and  r  and  their  partial  derivatives  are  equated  to  zero  to 
determine  unique  and  global  minimum  costs  in  each  case.   In 
particular,  they  claim  these  functions  are  strictly  convex 
in  the  case  where  lead  time  demand  (that  demand  occurring 


between  the  reorder  point  and  arrival  of  that  order)  is 
normally  distributed.   These  assertions  are  unproven 
however,  and  are  left  to  the  reader  as  a  series  of  exercises. 

The  first  exercise  asks  the  reader  to  show  that  a 
particular  terra  which  occurs  in  both  models,  say  J(Q,r) , 
is  convex.   It  is  obvious  that  the  other  terms  are  convex. 
It  would  then  follow  that  K(Q,r)  is  convex  in  Q  and  r  from 
additivity.   However,  it  was  pointed  out  in  1964  by  Veinott 
[Ref.  2]  that  J(Q,r)  is  not  always  convex.  In  1969,  Brooks  and  Lu 
[Ref.  3]  considered  the  convexity  in  the  backorders  model 
and  showed,  for  normally  distributed  lead  time  demand,  that 
J(Q,r)  is  convex  for  all  r    >_  \i ,    the  mean  lead  time  demand. 
In  1975,  the  nonconvexity  of  the  backorders  cost  function 
and  the  nonuniqueness  of  its  solution  was  demonstrated  by 
Minh  [Ref.  4]  for  the  case  where  lead  time  demand  is  normally 
distributed.   In  particular,  Minh  showed  through  counter- 
examples that  K(Q,r)  is  not  convex  and,  if  a  solution  exists, 
its  partial  derivatives  equated  to  zero  have  exactly  two 
solutions  vice  one  as  claimed  by  Hadley  and  Whitin.   Minh 
also  succeeded  in  classifying  one  of  the  solutions  as  a 
relative  minimum  and  the  other  as  a  saddle  point.   Further, 
Minh  showed  that  feasible  minimum  solutions  can  be  obtained 
where  r  <  y. 

As  Minh  has  suggested,  his  results  have  created  doubts 
as  to  Hadley  and  Whitin' s  claim  that  the  lost  sales  cost 
function  is  strictly  convex  for  the  case  of  normally  dis- 
tributed lead  time  demand.   Indeed,  the  only  mathematical 


difference  between  the  two  models  consists  of  an  additional 
term  in  the  inventory  carrying  cost  expression  of  the  lost 
sales  cost  function.   This  term  accounts  for  the  additional 
cost  incurred  by  carrying  as  inventory  that  portion  of  an 
arriving  lot  which  would  have  been  immediately  dispatched 
to  satisfy  accumulated  backorder  demand,  if  such  were  per- 
mitted.  This  is  consistent  with  the  authors'  assumption 
that  the  expected  number  of  lost  sales  is  equivalent  to  the 
expected  number  of  backorders,  all  other  factors  being 
equal . 

If,  as  will  be  shown  in  this  thesis,  K  is  not  convex 
as  claimed,  then  of  course  a  different  set  of  arguments  must 
be  provided  for  its  minimization.   In  particular,  the  con- 
straint imposed  by  assuming  a  nonnegative  reorder  point 
will  also  have  to  be  examined.   Hadley  and  Whitin  ignore 
this  constraint  in  their  development,  leaving  the  impression 
that  the  solution  is  always  to  be  found  in  the  region 
r  >  0,  Q  >  0.   Indeed,  they  are  not  clear  as  to  their 
intent  regarding  r  at  all.   In  the  beginning  of  their 
derivation,  they  require  r  to  be  positive;  in  the  exer- 
cises dealing  with  convexity  however,  r  is  allowed  to  be 
nonnegative.   Since  r  is  assumed  to  be  a  continuous  variable, 
there  seems  to  be  no  rational  basis  for  excluding  r  =  0  and 
it  will  be  assumed  herein  that  r  is  to  be  treated  as  a 
nonnegative  variable. 

It  is  the  purpose  of  this  thesis  to  clarify  the  issue  of 
convexity  for  the  lost  sales  case  and  completely  characterize 


its  solution,  including  the  effect  of  the  constraint  on 
r,  assuming  a  normally  distributed  lead  time  demand. 
Section  II  will  show  that  K  is  not  convex  and  that  the 
nonnegativity  constraint  on  r  precludes  Hadley  and  Whitin's 
claim  of  existence  in  that  the  simultaneous  solution  to 
the  partial  derivations  of  K  may  be  outside  its  domain  of 
definition.   Section  III  will  present  a  series  of  lemmas 
and  theorems  in  which  the  minimum  solution  is  completely 
characterized.   Section  IV  will  briefly  discuss  these 
results . 

For  a  normally  distributed  lead  time  demand,  Hadley  and 
Whitin  define  the  lost  sales  cost  function  as  follows: 


K(Q,r)  =  ^+TC[|+r-y]  +  [IC  +  TrA][(y-r)$(^+a<J>£=H.)] 


The  parameters  are  defined  as  follows: 

A  is  the  set-up,  or  ordering  cost. 

C  is  the  unit  cost. 

I  is  the  inventory  carrying  cost  per  unit  per  year. 

7T  is  the  cost  of  a  lost  sale. 

X  is  the  mean  annual  demand. 

a      is  the  standard  deviation  of  the  normally  distributed 
lead  time  random  variable. 

y   is  the  mean  of  that  random  variable. 

<p(z)   is  the  standard  normal  density  function. 

$(z)   is  its  complimentary  cumulative  distribution  function, 

00 

i.e.   ■?  ( z )  =  /  M  x )  dx . 
z 


For  ease  in  subsequent  algebraic  manipulation,  K(0,r) 
will  be  rewritten  as 


K(z,Q)  =  —  +  IC[£  +  az  +  n(z)]  + 


Q      l2    —        "^'J      Q 

r—  u 
where  z  =  — -  and  n(z)  =  a[<J>(z)  ~  z<f>(z)].   In  this  notation, 

the  set-up  cost,  inventory  carrying  cost,  and  lost  sales 
costs  are  clearly  identifiable.   The  domain  of  definition 
for  K,  as  a  function  of  z  and  Q,  is  defined  by  the  restric- 
tions z  >  -  ~,  Q  >  0. 
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II.   COUNTEREXAMPLES  TO  CONVEXITY 
AND  EXISTENCE  OF  SOLUTION 

A  function  of  two  or  more  variables  is  strictly  convex 
if  and  only  if  its  Hessian  matrix  is  everywhere  positive 
definite.   The  Hessian  matrix  is  positive  definite  if  and 
only  if  its  proper  diagonal  terms  and  determinant  are 
strictly  positive.   The  Hessian  for  K(z,Q)  is  derived  below 

The  Hessian  matrix  is  defined  as 


ft  = 


2 

3    K 

32K 

3z2 

3z3Q 

32K 

2 

3    K 

3Q3z 


3Q' 


Since  K(z,Q)  is  continuously  dif f erentiable  over  the 
region   z  >  -  —  ,  0  >  0 , 


2       2 

3^K     3  K 


3z3Q    3Q3z  * 


For  subsequent  derivations,  K  will  be  written  as 


K(z,Q)  =  ^[A  +  ttti(z)]  +  IC[|  +  az  +  n(z)]. 


Then, 


3K  _   TT.\ 

3z    Q 


n' (z)  +  ic[a  +  n*  (z) ] 
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where 


Thus, 


n ' (z)    =   -o$  (z)    . 

9k  ttAcj<J>(z) 


3z  Q 


+    ICa[l    -    *(z) ] . 


Hence, 


|f=   0  implies  t(>*)    =   QriIc°+  n  (1) 


J*(z)-  +   ica*(z 


3z2  2 


=  21^1  [QIC  +  ffX] 


2  2 

8    K  ^Aa$(z)  3    K 


3Q3z  2  3z30 


9K  A   r _  /    \  i  IC 

_    =     -    _[A    +     ^(z)]      +    __ 


Q 
Then, 


3K         n                   .  .                  _  .       .     /2A[A   +    irn  (z*)  ]  ,0, 

^  =    0  implies  0*    =\/ jg (2) 


92K         2ArA  ,    .  . 

[A    +    7TH  (z)  ] 


2  3 

3Q  Q 


The  Hessian  matrix  is  therefore 
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ft  = 


gtj)  (z)  [QIC  +  ttA] 


ttAg*  (z) 


L 


Q' 


it\o&  (z) 


2A  [A+  7TT1  (z)  ] 


Q' 


It  is  shown  in  the  Appendix  that  the  function  n(z) 
is  strictly  positive.   Thus  the  proper  diagonal  terms  of 
ft  are  strictly  positive  and  ft  is  positive  definite  if  and 
only  if  det  ft  (determinant  of  ft)  is  positive.   The 
determinant  of  ft  is  given  by 


det  ft  = 


7 


2cp  (z)  [A+  ttti  (z)  ]  [QIC  +  ttA]    -    tt2Aj$2(z) 


Define  the  expression  contained  within  the  outer  brackets 

as  V(z,Q).   Since  — j   >  0 ,  det  ft  >  0  if  and  only  if 

Q 

V(z,Q)  =  2  $  (z)  [A+  7m  (z)  ]  [QIC  +  ttA]  -  tt  Xa^>2(z)  >  0. 


With  this  result,  a  counterexample  to  the  claim  that  K(z,Q) 
is  strictly  convex  is  available.   The  parameters  utilized 
are  from  an  example  given  by  Hadley  and  Whitin  where  the 
lead  time  demand  is  normally  distributed  and  are  specified 
as  A  =  4000,  C  =  50,  I  =  0.20,  A  =  1600,  tt  =  2000,  a  =  50, 
M  =  750.   Then, for  the  policy  r  =  700  (z  =  -1.0)  and  Q  =  10, 
<Mz)  -  0.2420,  $(z)  =  0.8413  and  so, 
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V(z,Q)  =  2(.2420)   [4000 +  2000 (50) (.2420  +.8413)] 

x[(10)  (.20)  (50)  +  (2000)  (1600)] 
-  (2000)2  (1600)  (50)  (.8413)2 


=  1.7398214  x  1011  -  2.2649142  x  1011  <  0. 


Thus,  K(z,Q)  is  not  strictly  convex  throughout  the  region 
r  >  0 ,  Q  >  0  as  claimed. 

As  for  a  counterexample  to  the  existence  of  the  solu- 
tion, choose  the  following  parameters:   A  =  25,000, 
C  =  10,000,  I  =  0.99,  tt  =  10,  X    =  5,  a  =  10,  and  p  =  30. 
For  these  parameters,  the  iteration  procedure  described 
by  Hadley  and  Whitin  produces  a  solution  (z*,Q*)  =  (-3.10,5.06) 
It  may  be  verified  by  results  to  be  subsequently  derived 
that  V(z*,Q*)  is  positive.   It  will  be  observed,  however, 
that  the  reorder  point  r  corresponding  to  z*  =  -3.10  is 
r  =  z*a  +  u  =  -1.0  which  is  outside  the  domain  of  definition 
for  K(Q/r)  as  defined  by  Hadley  and  Whitin.   Basically, 
the  counterexample  is  a  consequence  of  the  fact  that  the 
solution  (z*,Q*)  is  independent  of  y,  the  mean  lead  time 
demand. 
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III.   CHARACTERIZATION  OF  THE  MINIMUM  SOLUTION;  PART  A 

It  is  convenient  to  separate  out  as  lemmas  certain 
isolated  results  that  support  the  main  results  presented 
in  this  thesis.   The  first  four  of  these  lemmas  are  stated 
in  this  section  and  derived  in  the  Appendix.   They  lead  to 
an  intermediate  result  stated  as  THEOREM  1,  a  theorem  that 
demonstrates  the  uniqueness  of  the  solution  (z*,0*)  to 
equations  (1)  and  (2)  and  establishes  that  solution  to  be 
an  unconstrained  relative  minimum  for  K. 

LEMMAs  5  thru  7 ,  also  stated  in  this  section  and  derived 
in  the  Appendix,  lead  to  THEOREM  2,  a  theorem  establishing 
that,  for  all  points  (z,Q)  where  z  f    z* ,    K(z*,Q*)  <  K(z,Q). 
The  main  result  will  then  be  presented  as  THEOREM  3  in  which 
the  minimum  solution  is  completely  characterized. 

Let  a  set  of  parameters  A,  C,  I,  it,  \,    a,    and  \i ,    all 
greater  than  zero,  be  given.   Define  the  following  functions 
for  -°°  <  z  <  °° : 


A(z)  .  ifMiuj .  ffT1(z) 

2IC[1  -  $(z)] 


G(z)     =    3(\)  (z)     +    z    -    z*(z) 


H(z)  ._*<*> 


[1  -  $(z) ]3 
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LEMMA  1:   Let  (z*,Q*)  be  a  solution  to  equations 
(1)  and  (2) .   Then, 

(a)   A  =  A(z*) 


^  J 


H(z 


ICQ 

ttA 


LEMMA_2:   (a)   The  function  G(z)  is  strictly  positive. 

(b)  The  function  H(z)  is  strictly  monotone 
decreasing. 

T*  C^  rr 

(c)  The  equation  H(z)  =  — *—  has  one  and  only 


ttA 


one  solution 


LEMMA_3:   The  function  A(z)  has  a  unique  minimum  at 

z   where  z   is  the  unique  solution  to  the 
o        o  ^ 


equation  H(z)  = 


ICj 
ttA 


LEMMA  4:   (a)   Lim  A(z)  =  0 


(b)   Lim   A(z)  = 

Z->— oo 


(c)   A(zq)  <  0 
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(d)   There  exists  one  and  only  one  z,  <  z 

1    o 

such  that  A(z,)  =  0.   Moreover,  if 


z  >  z,,  A(z)  <  0. 


THEOREM  1:   There  exists  one  and  only  one  (z*,0*)  which  is 

a  simultaneous  solution  to  equations  (1)  and  (2; 
and  that  solution  determines  a  relative  minimum 
for  K. 


PROOF : 


From  LEMMA  1,  let  z*  be  chosen  so  that 
A(z*)  =  A  >  0.   Then  z*  <  z,  from  LEMMA  4(d). 
From  the  strict  monotonicity  of  A(z)  for 
-co  <  2  <_  z,  ,  it  then  follows  that  z*  is  unique. 
Let  Q*  then  be  selected  according  to  either  of 
equations  (1)  or  (2).   Then,  from  LEMMA  1, 
(z*,Q*)  is  a  unique  solution  to  these  equations. 
To  show  that  K(z*,Q*)  is  a  relative  minimum, 
note  that  z*  <  z,  <  z  .   Then,  by  LEMMA  3, 
H(z*)  >  ^SrSL.   But  LEMMA  1(b)  establishes  this 
as  a  necessary  and  sufficient  condition  for 
positive  definiteness  of  the  Hessian  of  K  at 
(z*,Q*).   This  positive  definiteness  establishes 
that  K(z*,Q*)  is  a  relative  minimum  of  K  and 
uniqueness  has  already  been  established. 

Q.E.D. 
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III.   CHARACTERIZATION  OF  THE  MINIMUM  SOLUTION;  PART  B 

The  results  of  the  previous  lemmas  are  depicted 

graphically  in  Figure  1.   The  functions  A(z)  and  H(z)  and 

Ipq 
the  lines  A  and  -— -  are  shown  for  perspective  in  that  they 

represent  no  particular  set  of  parameters.   It  should  be 

noted  that  z*,  here  shown  negative,  may  be  positive  or 

negative. 


Figure  1 


It  has  been  established  in  PART  A  that  K(z*,Q*) 
constitutes  a  relative  minimum  and  that  (z*,0*)  is  unique 
In  this  part  it  will  be  demonstrated  that  K(z*,Q*)  does 


constitute  an  unconstrained  global  minimum.   However, 
it  was  shown  in  Section  II  that  the  constraint  r  >_  0 
(z  _>  -  — )  may  be  binding,  Therefore  that  constraint  must 
be  considered  in  the  minimization  of  K. 

For  the  succeeding  lemmas  and  theorems,  it  will  be 
convenient  to  define  the  following  functions  for  Q  >  0, 

-oo  <   Z   <  °°. 


Q]_(z)  = 


7T\$(Z) 


IC[1  "  $(Z) ] 


Q2(z)  =" 


2A  [A  +  irri  (z)  ] 


IC 


ZX(Q)  =  $ 


-1 


QIC 


QIC  +  ttX 


q^x(q: 


LEMMA  5 :   Let  z  be  fixed.   Then,  K(z,Q2(z))  <_   K(z,Q) 
for  all  Q. 


LEMMA  6:   If  z  <  z*,  Q2(z)  <    Q-i(z)  /  and  if  z  >  z*, 
Q2  (z)  >  Qx(z)  . 

LEMMA  7 :   Let  Q  >  0  be  fixed.   Then,  if  z  <  z-^Q)  , 
||  <  0  and,  if  z  >  zx(Q) ,  ||  >  0. 


THEOREM  2:   If  z  ^  z* ,    then  K(z*,Q*)<   K(z,Q)  for 
all  z  and  Q  >  0. 
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PROOF:       Let  z  be  fixed  and  v  =  (cos  9,  sin  8)  be  a 

unit  vector  from  (z,Q2(z))  in  the  direction  of 
(z*,Q*) .   The  directional  derivative  of  K  at 
(z,Q2(z))  may  be  written  as 


DK  v  -  i  cos  a  +  1  sin 


8K 

■*—  cos 

dZ 


9K 
since  -^  =  0  at  (z,Q2(z)).   The  proof  now 

proceeds  by  cases. 


Case  1:   z  <  z*.   Since  Q2(z)  i-s  strictly  monotone 
decreasing,  Q2(z)  >    Q2(z*)'  and  thus  -=-  <  9  <  2it. 
Hence,  cos  9  >  0.   Let  0,  =  0, (z)  and  Q?  =  Q~ (z) . 
By  LEMMA  8 ,  Q2  <  Q,  ,  in  which  case  z,(Q2)  >  z, (Q,) 
since  z,  is  strictly  monotone  decreasing. 
Since  z, (Q,)  =  z,  z, (Q2)  >  z  which  implies 

t-  <  0  by  LEMMA  9.   Hence,  D.,v  <  0  for 
z  <  z*. 


Case  2:   z  >  z*.   Since  0~(z)  is  strictly  monotone 
decreasing,  Q-(z)  <  Q_(z*),  and  thus  j   <    9  <  tt  . 
Hence,  cos  9  <  0.   Let  Q   =  Q,(z)  and  0   =  Q2(z). 
By  LEMMA  8,  Q2  >  Q,,  in  which  case  z^Q-)  <  z.,  (Q^) 

since  z.  is  strictly  monotone  decreasing.   Since 

3K 
z,  (Q,  )  =  z,  z,  (Q-)  <  z  which  implies  ^—  >  0  by 

11  1    Z  dZ 

LEMMA  9.   Hence,  D^v  <  0  for  z  >  z*. 
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In  all  cases,  D  v  <  0  so  that  K  is  strictly 
decreasing  from  any  point  (z,Q?(z))  toward 
(z*,Q*).   Thus  K(z*,Q*)  <  K(z,Q2(z))  and,  by 
LEMMA  5,  K(z,Q2(z))  <  K(z,Q)  for  any  Q  >  0. 

Q.E.D. 


THEOREM  3: 


If  z*  _>  -  -±,    then  K(z*,Q*)  is  the  global  minimum 


If  z*  <  -  j±,  then  K(-  |,Q2(-y/a)) 


value  of  K. 

is  the  global  minimum  value  of  K. 


PROOF: 


If  z*  >  -  — ,  then  the  corresponding  value  of  r 
is  nonnegative  and  r  is  in  the  domain  of  definition 
of  K.   From  THEOREM  2,  K(z*,Q*)  <  K(z,Q)  for  all 
z  ^  z*  and  Q  >  0.   Thus,  K(z*,Q*)  is  a  global 
minimum  if  Q  >  0  and  z*  >  -  — . 

—    (7 

If  z*  <  -  — ,  then  the  corresponding  value  of 
r  is  negative  and  is  thus  outside  the  domain  of 
definition  of  K.   Hence,  K(z*,Q*)  cannot  be  optimal 
However,  LEMMA  5  has  established  that,  for  any 
fixed  z,  K(z,Q2(z))  <  K(z,Q)  and,  by  THEOREM  2, 


K 


is  strictly  decreasing  from  any  point  (z,Q2(z 


toward  (z*,Q*)  .   Hence,  K(-  JpQ2(-  £)  )  1  K(z,Q2(z) 
for  any  z    -  —  (r    0) 


Thus ,  the  minimum  value 


of  K  occurs  on  the  boundary  z  =  -  —  and,  for 

z*  <  -  — ,  K  is  minimized  at  K(-  -,Q~.(-  — )  )  . 

a  o      z        o 


Q.E.D 
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IV.   DISCUSSION 

The  results  of  this  thesis  have  resolved  a  long-standing 
question  concerning  the  solution  to  the  lost  sales  case  of 
the  Hadley-Whitin  (Q,r)  model.   Since  the  model  is  widely 
used  by  practitioners,  it  is  an  issue  that  needs  resolution, 
for  the  entire  argument  given  by  the  authors  is  based  on  the 
supposed  convexity  of  the  objective  function.   That  not  being 
the  case,  it  is  necessary  to  supply  a  different  argument  for 
the  solution.   Convexity  aside,  the  boundary  constraint  can 
be  binding  and  appears  to  have  been  completely  ignored  by 
the  authors  in  their  development.   This  thesis  has  clarified 
both  these  issues. 

Beyond  the  mere  clarification  of  the  results,  this 
thesis  exemplifies  the  fact  that  attention  to  the  underlying 
mathematical  detail  and  rigor  cannot  be  ignored  in  modeling 
a  problem.   Fortunately,  in  this  case,  the  solution  remains 
partially  valid  but  for  a  different  set  of  reasons.   This 
could  only  have  been  discovered  by  a  close  and  detailed 
examination  of  the  mathematics  of  the  model.   In  the  process 
of  that  examination,  the  role  of  the  boundary  condition  has 
revealed  itself  and  allows  for  complete  characterization  of 
the  solution. 

As  for  recommendations  for  further  study,  it  should  be 
pointed  out  that  only  the  case  of  normally  distributed  lead 
time  demand  has  been  considered.   There  are,  however,  many 
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other  assumptions  that  might  be  made  concerning  this  demand. 
Among  those  of  special  interest  to  Navy  applications  are 
those  of  the  Poisson  and  negative  binomial  distributions. 
An  investigation  similar  to  that  presented  here  for  at  least 
these  two  cases  would  seem  to  be  called  for  and  is  a 
recommendation  of  this  thesis. 
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APPENDIX  A 
PROOFS  OF  LEMMAS 

This  Appendix  contains  the  proofs  of  LEMMAs  1  thru  7 
inclusive.   Prior  to  their  proofs,  it  will  be  convenient 
to  prove  four  auxiliary  lemmas  which  are  appealed  to  in 
the  main  argument.   For  clarity,  the  auxiliary  lemmas  will 
be  numbered  with  the  letter  "A"  following  the  number. 

LEMMA  1A :   The  function  n(z)  is  strictly  monotone  decreasing 
for  -oo  <  z  <  °°  and  everywhere  strictly  positive 
over  that  region. 

PROOF:      n  (z)  =  cru(z) 

where 

u(z)  =  $  (z)  -  z$ (z)  . 

Now, 

u'  (z)  =  -z<j>(z)  -  z(-<f>(z)  )  -  *(z) 

=  -  $(z)  <  0      for  all  z. 
u"  (z)  =  4>  (z)  >  0       for  all  z. 
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Clearly, 


Lim      u(z)    =    Lim       (cj)  (z)  -z<S>  (z)  )    =    0    +   °°  =    +°° 

2->— oo  Z-*—  °° 


And 


Lim  u(z)    =    Lim    («|>(z)-z$(z))    =    -Lim   z$(z) 

2  -*00  Z  "*-0°  Z  ~>"oc> 


=    -    Lim 


2->oo 


$(z) 


Since  this  limit  results  in  the  indeterminant 

oo  ^ 

form  — ,  L'Hopital's  Rule  is  applied  repeatedly 


Lim  — i =  -  Lim 


2  -h»   2  ~>0° 

$  (z) 


1 

T"               ^2(Z) 

-    Lim    ? j  \ 
d>  (z) 

4>  (z) 

*2(z) 

,.   -2(|>(z)3(z) 

Lim  —7 — c 

-z4>  (z) 
2-^-00 


_  .   2$(z)  _  n 

Lim  -  0 . 

z 

z->-co 


Thus  n(z)  =  au(z)  is  a  strictly  convex  decreasing 
function  approaching  0  as  z-h-°°  and  hence 
H (z)  >  0  for  all  z . 
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^^Jh--     QxW  =  icTT^fW 


is  strictly  monotone 


decreasing  and  everywhere  strictly  positive 


PROOF 


Bi(«)  -& 


-$(z)(Hz) 

[1  -  $(2)  ] 


+ 


-<j>(z) 


2     [1  -  $(2)  ] 


—  irXcf)  (z) 


IC[1  -  $(2)  ] 


<  0     for  all  z. 


Q,  (2)  is  everywhere  strictly  positive  by  inspection 


Q.E.D 


LEMMA  3A:   Q2<z)  = 


2A  [A  +  ttti  (z)  ] 


IC 


is  strictly  monotone 


decreasing  and  everywhere  strictly  positive 


PROOF:      Squaring  Q~(z)  and  differentiating  implicitly 


QJ(z)  = 


\ttg$  (2 ) 
Q2(z)IC 


<  0     for  all  z 


Q-(z)  is  everywhere  strictly  positive  by  inspection 


LEMMA  4A:   z.,(Q) 


=  * 


-1 


QIC 


QIC  +  ttA 


Q.E.D. 


is  strictly  monotone 


decreasing. 
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PROOF 


Qx(z) 


7TA.$  (Z) 


IC[1    -    $(z)] 


Thus, 


(z)    = 


OIC 


QIC    +    ttA 


or, 


LEMMA    1 


z    = 


,-1 


QIC 


QIC     +     TTA 


=    Zl(Q)     =    Q~1(Q) 


Since  z, (Q)  is  the  functional  inverse  of  Q, (z 
and  Q, (z)  is  strictly  monotone  decreasing, 
z, (Q)  is  strictly  monotone  decreasing. 

Q.E.D. 


(a)   A  =  A(z*) . 

\=  \      0     if  and  only  if 


H(z* 


iCg 


Proof  (a) .    Equating  the  two  expressions  for  Q*  given  by 
equations  (1)  and  (2)  yields 


IC[1  -  l>(z*)  ] 


2\ [a  +  ^n  (z*) ] 
IC 
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or, 


2       2 
/^    (z*) 


IC[1    -    $(z*) ] 


J       =        2  [A    +     TTTl(Z*)  ] 


Solving    for   A, 


7T2A$2  (Z*) 


2IC[1    -    $(z*) ] 


-     1TT1  (Z*)         =        A 


that   is, 


A(z*)    =   A 


Proof     (b)  .       V(z*,Q*)    =    2cf)  (z*)  [A+ttti  (z*)  ]  [Q*IC+ttA]    -    TT2\a$2(z*) 


Substituting  A(z*)  for  A  from  part  (a)  and 


ttA©(z*) 


Q*   =   -, 


IC[1    -    $(z*) ] 


from  equation   1, 


V(z*,Q*)    =    2((>(z*) 


2       2 

TT      A<£>      (  Z* 


2IC[1    -    $(z*) ] 


ttA<I>(z*) 


[1    -    $(z*) ] 


TT  A 


2         2 
-    it    Aa$    (z*) 


=    2<J>(z*) 


2       2 
T    A$     (z*) 


2IC[1    -    $(z*) ] 


ttX 


[1    -    $(z*)] 


2         2 
-    TT    Aa$    (z*) 
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V(z*,Q*)  =  tt2A$2(z*)    JLhlLlll -  0 

IC[1  -  $(z*) ] 


Since  it  X*  (z*)  >  0,  V(z*,Q*)         0   if  and 


only  if 


-  a      {  =  > 


IC[1  -  #(Z) ]3 


ic 

Since  -y  >  0,  this  condition  may  be  written  as 


=  H(z*)  lCa 


[1  -  *(z*)]3 


Q.E.D. 

LEMMA  2 :   (a)   The  function  G(z)  =  3<(>(z)  +  z  -  z$(z)  is 

strictly  positive. 

(b)  The  function  H(z)  =  — ^-^ *■  is  strictly 

[1  -  *(z)] 

monotone  decreasing. 

rco 

(c)  The  equation  H(z)  =  — r—  has  one  and  only 

TT  A 

one  solution. 

Proof  (a).   The  proof  proceeds  by  cases. 

Case  1.   z  _>  0. 

G(z)  =  3<J)(z)  +  z  -  z$(z)  =  u(z)  +  2$(z)  +  z 
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Since  u(z)  =  <|>(z)  -  z$(z)  is  strictly  positive 
for  -co  <  z  <  «,  G(z)  >  0. 

Case  2.   z  <  0.   Let  z  =  -i"  where  z~  >  0 . 
Then, 

G(z)  =  G(-z)  =  3  0  (-z)  -  z  +  z~$(-z) 

=  3<j>(z)  -  z  +  z[l  -  *(z)  ] 

=  3<j>(z")  -  z$(z) 

=  20  (z)  +  u(z)  >  0. 

Proof  (b) .   Differentiating  H(z), 


2 

Hi  (z)  _  -3  0  (z) z<Hz) 

[1  -  $(z) ]4    [1  -  $(z) ]3 


30(z)    , 


[1  -  $(z) ]3   [1  -  $(z)  ] 


30(z)  +  z  -  z$(z 


[1  -  $(z) ]4 


The  expression  in  brackets  is  G(z)  which  has  been 
shown  in  part  (a)  to  be  strictly  positive.   Hence, 
H' (z)  <  0  and  H(z)  is  strictly  monotone  decreasing 


Proof  (c) .   Since  Lim  $(z)  =  0, 
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Lim  H(z)  =  Lim  — li^J =  o 

Z^oo  z->oo    [i  _  $  (Z)  ] 


Lim   H(z)  =  Lim   H(-z)  =  Lim 


3 

z ->■-<»        z->-°°  z-*-°°   $  (z) 


This  leads  to  the  indeterminate  form  ^r.   Applying 


L'Hopital's  Rule, 


Lim  4ill_  =  Lim     -*♦<«> 


3         "^au     2 
Z->co   $  (z)      2->°°   -3$  (z)cj>(z) 


Lim 


2 

2-voo    3$   (z) 


Since  H(z)  is  everywhere  positive  and  strictly 
monotone  and  — =r—  >  0,  the  equation  H(z)  =  — «— 

77  A  ^  77  A 

can  have  one  and  only  one  solution. 

Q.E.D. 


LEMMA  3:   The  function  A(z)  has  a  global  minimum  at  z 
^  o 

where  z   is  the  solution  to  the  equation 
Proof:     It  is  convenient  to  define  the  following: 


,  .      $  (z)  ttA 

g(z)  =  ■ ,      a  =  —  • 

[1  -  $(z)] 
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A(z)  can  then  be  written  as 


A  (Z)      =     7T 


a      2  , 


)    -  n(z) 


Then, 


A'  (z)    =  7T 


a   g(z)    g1 (z)    +    a$ (z) 


Since   g1 (z)    = 


-<j>(z) 


[1   -    $(z)] 


A*    (Z)      =     7T 


-ttA$(z)  4>(z) 
IC[1    -    *(z)] 


+   a$ (z) 


7T    \$(z) 
IC 


ICj  0  (z) 

ttA 


[1    -    <D(z)] 


Since  *  i,    >  0,  A"  (z)    {  =  ^    0   precisely  when 


H(z) 


ICa 


By  LEMMA  2,  there  exists  a  unique  z   such  that 


ho  )  =  £& 

o     ttA 
of  H(z) , 


Moreover,  from  the  strict  mono tonicity 


H(z) 


ICa 

uA 


for 
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It  then  follows  that 


{=)        0     for     z    l=\        z 


establishes  z  ,  the  only  critical  value,  as 
yielding  a  unique  minimum  for  A(z) . 

Q.E.D. 


LEMMA  4:   (a)   Lim  A(z)  =  0 


(b)   Lim   A(z)  = 

2->— co 


(c)   A(z  )  <  0 
o 


(d)   There  exists  one  and  only  one  z,  <  z 

i    o 

such  that  A(z,)  =  0.   Moreover,  if  z  >  z, 
A(z)  <  0. 


2  2 

ir      X    $    (z) 
Proof    (a).      A(z)    =  J — j  -    7rau(z 

2IC[1  -  $(z)  ] 


It  has  already  been  shown  that  Lim  u(z)  =  0 

z-*-<» 

Hence, 


2   2 

Lim  A(z)  =  Lim  * — • j  ~    Llm  ^cufz) 

z^00        z-*°°2lC[l-$(z)]     z^-00 


=  0. 
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Proof    (b) .       A(-z)    =   — 


Tra 
a 


(-z) 


[1    -    $(-z)] 


-    a  u(-z) 


where 


a   = 


2ICq 


A(-z)    =    ^ 


■no 

a 


[1  -   Hz)  ] 
$2(z) 


-   a  [<fr  (-z)    +    z$(-z)  ] 


a 


r  2 

[i  -  $(z)r 

$2(z) 


-   a  [u  (z)    +    z] 


a 


[1    -    $(z)] 
$2(z) 


-  az 


-   TTcru(z) 


Since   Lim     A(z)    =    Lim  A(-z)       and      Lim  u(z)    =    0, 


z-*-°° 


Lim     A(z)    =   Lim 


7TQ 


2->oo 


[1   -    Hz)  ] 
$2(z) 


-    az 


TTO        T  ■  [1 

—      Lim      - — 
a 


2  2 

$ (z) ]       -    az$    (z) 

*2(z) 


^2  ,    . 
But,    Lim   z$     (z)    = 


Lim 


1 


$      (Z) 


=      Lim 

2->-co 


2$(z) 
$3(z) 


1  t^   Lili 
"t"  Lim  — r-7 — r 

2  4)  (z) 
z-*-<» 


1  T.  -3^(z)(Hz) 
7T  Lim  r~7 — r 

2  -z4>  (z) 
z-*-°° 


3      Tirr,  *1L*1 

■x-  Lim  

2  z 

Z->co 


=    0 
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2 
and   Lim    [1   -    $(z)]      =   1.      Hence, 


2-H30 


Lim     A(z)    =    Lim  A(-z)    =   °°, 

2->--oo  Z-^00 


Proof  (c) .   The  proof  is  by  contradiction.   Assume 

A(z  )  _>  0.   Since  A(z)  is  strictly  increasing 

on  the  interval  z  >  z  ,  A(z)  >A(z)  >0  for  all 

o  o      — 

z    >    z    .       This    implies    Lim  A(z)     >    0    contradicting 
o 

z->°° 

the  fact  that  Lim  A(z)  =  0  established  in  part  (a) 

2->oo 

Proof  (d) .   It  has  been  shown  in  LEMMA  3  that  A' (z)  <  0 

for  -°°  <  z  <  z   and  in  part  (b)  above  that 

Lim  A(z)  =  °°.   Since  A(z  )  <  0,  there  must  exist 

o 

2->—  oo 

-  oo  <  2,   <  z   such  that  A(z,)  =  0.   That  z,  is 
1     o  1  1 

unique  in  the  interval  -OT  <  z    z   follows  from 
the  strict  monotonicity  of  A(z)  in  this  interval. 

Q.E.D. 


LEMMA  5 :   Let  z  be  fixed.   Then,  K (z ,Q2 (z ) )  <  K(z,Q)  for 

all  Q. 


Proof:     E'or  any  fixed  z, 


3K  _  -  a[A  +  ttti  (z)  ]    IC 


oQ  r2  2      ' 


II  =  0     implies     Q  ^\J ~ ~ fc =  Q2(z) 
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9  K    2X[A  +  Trn(z)]  .  n 
— j , >  0. 


3Q 


Q' 


Thus,  K(z,-)  is  a  convex  function  with  a  unique 
minimum  at  Q~ (z) . 


Q.E.D. 


LEMMA  6:   If  z  <  z*,  Q2(z)  <  Q,(z)  and  if  z  >  z*,  Q2(z)  >  Q,(z) 


Proof:     Define  the  following  functions  for  -«  <  z  <  <»: 


D(z)    =    W2(z)     -    W1(z) 


where 


Wx(z) 


Ql.(z) 


IC 


2  2 

$    (z) 


[1   -   <Mz)  ] 


and 


W2(z) 


Q^(z) 


2  A  [A    +     TTTl  (z)  ] 
IC 


Then, 


D' (z)    =   W2(z)    -    W^(z 


-2tt  Agfc  (z 
IC 


+    2 


ic 


$(z)    H(z) 


=    2 


ic 


<S>(z) 


H(z)     - 


ICa 
ttA 
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Thus,  by  LEMMA  2, 

V 

D' (z)    \=)         0     if  and  only  if 


z  . 
>■     ° 


Hence,  z   is  a  relative  maximum  for  D.   Observe 
o 

that,  since  Q,(z)  >  0  and  Q2(z)  >  °# 


Q2(z)    M    Q1(z)    iff   W2(z)    <=j>   W1(z) 
iff   D(z)    \=}         0    and  D  (z)  =  0  if  and  only 


if  z  =  z*. 

Suppose  z  <  z*.  Since  z*  <  z   and  D  is 

increasing  on  (-<*>, z  ),  D(z)  <  D(z*)  =  0.   Hence, 
Q2  (z)  <  Q1(z)  . 

Suppose  z  >  z*.  If  z*  <  z    z  ,  then 

0  =  D(z*)  <  D(z)  or  Q2(z)  >  Qi(z)-   If  z   <  z, 
then,  since  D  is  decreasing  on  (z  ,°°)  ,  and 

lim  D(z)  =  ^^  >  0,  D(z)  >  0   or   Q2(z)  >  Q1(z). 

Q.E.D. 

LEMMA  7:     Let  Q  >  0  be  fixed.  If  z  <  z, (Q) ,  then  4^  <  ° > 

L  a  Z 

while  if  z  >  z, (Q) ,  |£  >  0. 

I  dZ 

Proof:       For  any  fixed  Q, 


=  iCa  -  [IC  +  — ]a$(z 

dZ  Q 
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3K 

3z 


=  0 


implies 


z  = 


-1 


QIC 

QIC  +  irX 


=   z2(Q) 


2 

77 


=  [IC  +  TTA]a(t>  (z)  >  0 


Thus,  K(*,Q)  is  convex  with  a  unique  minimum 

9K 
at  z,  (Q)  .   Hence,  if  z  <  z,  (Q)  ,  ^r—  <  0,  while 

_L  J.        d  Z 

if  z  >  z1(Q),  ||  >  0. 


Q.E.D. 
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